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Abstract: We develop the "generalized Gordon ansatz" for the ghost-free versions of 
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to greatly simplify calculations. The ansatz allows explicit calculation of the matrix 
square root 7 = \/g~ l f appearing as a central feature of the ghost-free analysis. In 
particular, this ansatz automatically allows us to write the effective stress-energy tensor 
as that corresponding to a perfect fluid. A qualitatively similar "generalized Kerr- 
Schild ansatz" can also be easily considered, now leading to an effective stress-energy 
tensor that corresponds to a null fluid. Cosmological implications are considered, as 
are consequences for black hole physics. Finally we have a few words to say concerning 
the null energy condition in the framework provided by these ansatze. 
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1 Introduction 



The ghost-free massive and bimetric gravities have recently attracted considerable at- 
tention. Early (pre-ghost-free) work included [1-3], while notable contributions to the 
recent ghost-free discussion include [4-12]. Recent work on the relation between mas- 
sive gravity and bimetric gravity has been reported in [13], and the status of the null 
energy condition (NEC) is explored in [14]. Cosmological implications are discussed 
in [15-21], while the status of black holes is considered in [22-29]. 

An overarching theme of all the ghost-free technical calculations is the need for 
manipulating and working with the matrix square root 7 = \Jg~ l f built out of the 
foreground and background metrics. Some particularly useful mathematical results 
related to dealing with this quantity are reported in [14]. In the current article we 
shall investigate a specific ansatz for the relation between foreground and background 
metrics that is sufficiently general to be able to represent almost all spacetimes of 
physical interest, while being simultaneously sufficiently special to make computations 
simple. The ansatz is inspired by the "analogue spacetime" programme [30-32], but 
has wide applicability and validity outside of that original context. (Note that the 
analogue spacetime programme could play a relevant role in this context not only as an 
inspiration, but even as a scenario where bimetric models could naturally emerge [33].) 

The paper is organized as follows: The general framework for ghost-free massive and 
bimetric gravity is very briefly summarized in section 2. Then in section 3 we introduce 
a particularly powerful ansatz which can be obtained by generalizing the Gordon metric 
(normally arising in relativistic optics). This "generalized Gordon ansatz" allows us 
to write the effective stress-energy tensor as that of a perfect fluid, see subsection 3.2. 
Despite the rather strong conclusions implied by this ansatz, it is still general enough to 
describe the most common physically interesting spacetimes considered in the literature 
- see in particular subsection 4. In section 5 we present the "generalized Kerr-Schild 
ansatz" and discuss the implied stress-energy tensor, now that of a null fluid. We 
conclude with a discussion in section 6. 

2 Stress energy tensor in bimetric gravity 

The action of bimetric gravity can be expressed quite generally as [13, 14] 



S 



1 



d A x^g~ {R(g) + 2 A - 2m 2 L int (g, /)} + S (m) 



lQnG 




(2.1) 
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Here S( m ) an d S'(m) are the usual matter actions, with the foreground and background 
matter fields being coupled only to the foreground and background metrics and 
f^u, respectively. All interactions between these two sectors are confined to the term 
Lint(g, f), which is an algebraic function of g and /. We can recover the action of mas- 
sive gravity by considering k — e — [13]. In this case, we would have an aether theory 
in which the dynamics of the physical metric g^ v depends on a now non-dynamical 
background metric f^ u . (Up to this point, the analysis is identical to the pre-ghost-free 
analyses of the 1970s [1, 2], and later [3].) 

As is now well-known [4], the action (2.1) is ghost-free if and only if the interaction 
term can be written as a linear combination of the elementary symmetric polynomials 



of the eigenvalues of the matrix 7 = \J g 1 f. In view of the results of [14] we can write 

Uu = l c \ 9a P Yu] that is / = 7 T g 7. (2.2) 
We can then, in 3+1 dimensions, express the interaction Lagrangian as [13] 

Lim = 011 ei(7) + "2 e 2 (j) + a 3 e 3 (~/), (2.3) 

where as usual 

ei(7) = tr[ 7 ]; (2.4) 

e 2 ( 7 ) = ^(tr[ 7 ] 2 -tr[ 7 2 ]); (2.5) 

e 3 ( 7 ) = ^(tr[ 7 ] 3 -3tr[ 7 ]tr[7 2 ] + 2tr[7 3 ]). (2.6) 

The two remaining non- vanishing polynomials, eo( 7 ) = 1 and 64(7) = det( 7 ), have 
been absorbed into the foreground and background cosmological constants. The three 
parameters appearing in the interaction term are not fully independent. If one addi- 
tionally requires that L int should (for weak fields) be of the canonical Fierz-Pauli form, 
then [7, 13] 

at + 2a 2 + a 3 = -1. (2.7) 

Both metrics have exactly the same status. The interaction term satisfies the reciprocity 
relation [4, 13] 

4 4 
y/=g Lint (7) = \ fZI 9 ^2 a i e i(l) = V~f ^2 a ^-i ei(7~ X ) = L int (7~ 1 ). (2.8) 

8=0 1=0 

That is, the entire theory could be equivalently re-expressed using / as the foreground 
metric and g as the background. Varying the action (2.1) with respect the two metrics, 
we obtain two sets of equations of motion [13]: 

G% - A 6» u = m 2 T» u + 8ttG T<- m K, (2.9) 
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and 

k (G* v - A 5"„) = m 2 T^j/ + e 87rG T (m)/ V (2.10) 

Here 

T^ = r^-5^L int , (2.11) 

and 



with 



T"„ = -^=| t"„, (2.12) 

^ = 7%^. (2.13) 

The indices of equation (2.9) and (2.10) must be raised and lowered using g and /, 
respectively. The equations of motion of the g-space (/-space) are modified with respect 
to those of general relativity by the introduction of an effective stress- energy tensor 
associated to the interaction between the two geometries (and the quantity re/e). It is 
a slightly non-trivial exercise, see [14], to verify that the effective stress-energy tensors 
are indeed symmetric. In massive gravity with a non-dynamical background one would 
only have the first set of equation of motions (2.9), see [13, 19]. Both effective stress- 
energy tensors fulfill the Bianchi-inspired constraints 

V M T^ = 0; V„T M „ = 0. (2.14) 

Once one constraint is enforced, for example V^T^V = 0, then the other is also auto- 
matically fulfilled [19, 20]. 

The modifications to the two sets of equations of motion are very closely related. 
Everything can be expressed in terms of a single mixed-index tensor r^, which can be 
written as a cubic polynomial in the matrix 7. Specifically [14] 

t» v = (on + «2 e a (7) + a 3 e 2 ( 7 )) -f u - (a 2 + a 3 61(7)) { 7 Y, + a 3 { 7 3 } M ,- (2.15) 

This yields the effective stress-energy tensor as viewed by the g-space, up to a cosmo- 
logical constant contribution coming from second term in the RHS of equations (2.11). 
Explicitly 

r^ = -(a 1 e 1 (7) + « 2 e 2 (7) + « 3 e 3 (7)) 8\ + t» v . (2.16) 
We also note the foreground-background symmetry 

yf^jT^ + y/^fT,, = -V=gL iat (j) 8\ = -7=7 Wt -1 ) (2.17) 
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3 Generalized Gordon ansatz 



Let us now consider the following ansatz relating the two metrics: 

/ F ^ 2 {V + (W (3.1) 

Here V 1 is a unit timelike vector with respect to g^, so that g l _ lv V^V u = — 1, and 
we set = g^vV 1 ' ■ Furthermore Q and £ are arbitrary functions. Noting that both 
metrics must have Lorentzian signature, we see £ < 1. So we can write £ = 1 — ( 2 , 
with C 7^ 0. Furthermore, without any loss of generality we can choose ( to be strictly 
positive, that is ( > 0. Then 

= n 2 {( 9flu + vm - c 2 v,m,}. (3.2) 

Mathematically, this ansatz can be thought of a conformal transformation Q, combined 
with a stretch ( along the timelike direction parallel to V. Note that (see figure 1) the 




Figure 1. Light cones for the generalized Gordon ansatz: Depending on the value of the 
parameter £ the foreground light cones lie strictly inside, on top of, or strictly outside the 
background light cones. 

relative position of the light cones is very strongly correlated. 
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• For < ( < 1 the light cones of / lie strictly inside the light cones of g. 

• For ( = 1 the light cones of / and g coincide; the two metrics are conformally 
related. 

• For C > 1 the light cones of / lie strictly outside the light cones of g. 
3.1 Historical background 

Physically, this ansatz is a significant generalization of Gordon's 1923 metric, which was 
originally developed within the context of geometric optics [34]; see particularly [30, 31] 
for recent discussions. The original Gordon metric, which describes ray optics in an 
optical medium of 4-velocity and constant refractive index n in Minkowski space, 
is expressed as 

G»u = (Tfox, + V/jYv) - ^ V/jYu, (3.3) 

where = r\ ilv V v . Our ansatz (3.2) corresponds to a generalization of (3.3). One 
should: 

(i) Generalize the refractive index n to be a position-dependent function (1/n 2 — > ( 2 ); 

(ii) Introduce a position-dependent conformal factor Q, (which does not affect the 
light cone structure); 

(iii) Consider a non-trivial input metric (g^ u instead of r]^); 

(iv) Relabel the output metric (f^ instead of Q^). 

For this reason, we call our ansatz the "generalized Gordon ansatz". The ansatz can 
also naturally be inverted to express g^ u in terms of f^: 

g, u = Q- 2 {(f, u + V,V U ) - C 2 V,V U }, (3.4) 

where V is now normalized in terms of f] lu V fM V V = —1, and we set V ^ = f^yV"- For 
consistency we must then enforce V ^ = an d V = Vl^ 1 ^ 1 V 1 . 

This generalized Gordon ansatz also has exactly the same structure as the (geomet- 
ric) relativistic acoustic metric of [35-37]; see particularly [30, 31] for recent discussions. 
Our ansatz, as expressed in equation (3.2), exactly agrees with the form of the geomet- 
ric acoustic metric appearing in reference [35-37]. (Here we have denoted the second 
metric f^ u instead of Q^ v . Furthermore in the present situation ( 2 is not necessarily 
restricted to be less than one. In contrast, in [35, 37] ( 2 = c 2 /c 2 < 1, simply because 
the speed of sound is subluminal, though "tachyacoustic cosmologies", based on [36], 
are now fashionable.) 
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3.2 Perfect fluid effective stress-energy tensor 

Using this ansatz we can express the matrix 7 in a simple form by noting that 

iTU = V 2 {{5\ + VV V ) - ( 2 VV V } = ( 7 2 )V (3.5) 
Then the matrix square root is easy to extract 

= (y/fF^Yv = + V^V V ) - C V>V V }. (3.6) 

Furthermore 

( 7 n )"„ = Q n {(8^ + V»V U ) - c V»V V }. (3.7) 

The symmetric polynomials appearing in the interaction term (2.3), given by equations 
(2.4), also take on a particularly simple form. They are 

ei ( 7 ) = Q{3 + O; e 2 ( 7 ) = 3Q 2 {1 + O; e 3 ( 7 ) = ^ 3 {1 + 3C}. (3.8) 

Note that these expressions are linear in (. This can be traced back to the fact that 
only one of the eigenvalues of 7 is at all ^-dependent, and that in a linear fashion. 
Taking into account expressions (3.5), (3.6), and (3.8), we see that equations (2.15) 
and (2.16), imply that the stress-energy tensor takes a perfect fluid form: 

T» v = V {8\ + V^Vy) + p V»V V . (3.9) 

Here we have defined 

p = tt (3«a + 3a 2 Vt + a 3 tt 2 ) , (3. 10) 

and 

p=-n [(2ai + a 2 n) + {oci + 2a 2 tt + a 3 n 2 )(] • (3.11) 

Note that these expressions are also linear in (, ultimately for the same reason as above. 
It is also useful to write 

T%=p(5^ + V^V u )+pV^V v , (3.12) 

defining 

p = -Q( («i + 3a 2 tt + 3a 3 tt 2 ) , (3.13) 

and 

p = n [ai + a 2 n(2 + C) + a 3 n 2 (l + 20] . (3.14) 

Here we have used the fact that T fJi u and r M „ only differ in a cosmological constant 
contribution. Moreover, we have 

p+p=p+p= ft( ai + 2a 2 Vl + a 3 n 2 )(l - C). (3.15) 
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Note that ( = 1 (£ = 0) corresponds to a cosmological constant, as should be expected, 
since in this case ansatz (3.1) reduces to a conformal relation between the two metrics. 

Moreover, due to the Bianchi-inspired constraint (2.14), the stress-energy tensor not 
only takes the perfect fluid form (3.9), but also obeys the usual conservation equations. 
These are 

V"*V M p+(p + p)V>% = 0, (3.16) 

and 

(p + p)V%K + ( gfiV + ^K)V> = 0. (3.17) 

On the other hand, in view of equations (3.10) and (3.15), it can be noted that the 
effective fluid automatically fulfills the additional constraint: 

P + P V^Q 

V,p = 3 — — . (3.18) 

From equations (3.16) and (3.18) one can conclude that, once Q and V a have been 
fixed in (3.2), the following equation is automatically fulfilled 

Finally, the effective stress-energy tensor in the /-space (2.12) also takes on a perfect 
fluid form 

T\ = p (<5"„ + V^V U ) +pV^ V v . (3.20) 

Now V ^ = Q( Vp, as defined in the previous section, is normalized in terms of 
fiivV^V 1 ' = — 1, and so 

V = -^k P] ~ P = ~^k P - (3 - 21) 

(This last step could equivalently be obtained by noting \J—f = fi 4 C \f~9-) Perhaps 
more tellingly 

P + P = -^ (P + P) = -^ (P + P)- (3-22) 

Thus, the sum of the effective energy density plus the pressure of the effective perfect 
fluid due to the modification of general relativity in the background space has opposite 
sign to that characterizing the effective perfect fluid in the foreground space. (This is a 
special case of the more general NEC- violating behaviour discussed in [14] . For general 
background on the NEC and its variants, see [38-41].) Some of the possible conse- 
quences of this sign flip are discussed in the particular case of cosmological solutions in 
section 4. 
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4 Generality of the ansatz 



Up to now, we have shown that the ansatz expressed through (3.1) or, equivalently, 
through (3.2), implies that the bimetric-induced effective stress-energy tensor takes the 
form of a perfect fluid. Moreover, the ansatz constrains the effects of the modifica- 
tions of the theory in the two spaces so tightly that (apart from the "trivial" case of 
pure cosmological constant) the violation of the NEC condition in one space follows 
straightforwardly from the fulfillment of the NEC in the other space [14] . Nevertheless, 
the ansatz is general enough to describe many different physically relevant situations. 
Some particularly relevant examples are considered below. 

4.1 Cosmological solutions 

Let us now assume that the metric takes the form of a FLRW spacetime in a certain 
coordinate patch. This is 



dsi 



-dt 2 + a(tf 



dr 2 



1 — kr 2 

Application of the ansatz (3.2) leads to 

dr 2 



+ r 2 {d9 2 + sin 2 9 dcf) 2 } 



(4.1) 



dsj = -n(ty U(t) 2 dt 2 + a(ty 

This metric can be written as 

ds) = -N(t) 2 dt 2 + b(t) 2 



kr 2 



dr 2 



1 — kr 2 

where the scale factor and lapse function of metric / are 

b{t) = Q(t)a(t), 

and 



+ r 2 {d9 2 + sin 2 fl d(f> 2 } 



+ r 2 {d6 2 + sin 2 9 d(f> 2 } 



(4.2) 

(4.3) 

(4.4) 
(4.5) 



N(t) = n(t)t(t), 

respectively. Thus, the assumption of the generalized Gordon ansatz (3.2), when con- 
sidering a cosmological foreground metric, implies a cosmological background metric 
with a different scale factor and cosmic time, as should be expected. The only re- 
strictions coming from forcing that relation is that the two metrics should both be 
diagonal in the same coordinate patch, and have the same sign of spatial curvature 
k, conditions that certainly should be expected to be fulfilled by general cosmological 
solutions [20, 21]. 
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The conservation equations (3.16) take on the usual cosmological form. Thus 

P + S(p + p) a =0. (4.6) 
a 

Taking into account equation (3.15) and the derivative of equation (3.10), this implies 

§ = -(l-0 d . (4.7) 

So, defining the Hubble parameters for both metrics in the usual way, H g = a /a and 
Hf = b/b, and taking into account the definition of b (4.4), we see that the two Hubble 
parameters are very tightly intertwined 

H f = C H g . (4.8) 

This is equivalent to the relation obtained in reference [21] to characterize their "branch 
two" solutions, the general class of cosmological solutions in bimetric gravity. 
Moreover, rewriting (4.7) as 

Q a .a j . 

o + ~ = C- (4-9) 



in view of (4.4) we see 



whence 



l = ( A , (4.10) 
b a 

b a dlnb/dt . 

( = = — . (4.11) 

a b dlna/dt 

Taking into account (4.5) and (4.4), this implies a lapse function 

N(t) = -. (4.12) 
d 

So, in this case, one can write the metric of the background space (4.3) using only one 
arbitrary function b(t): 

2 



dsi = - I t I df + b(tf 



a 



dr 



2 



1 — kr 2 



+ r 2 {d6 2 + sin 2 9 dcj) 2 } 



(4.13) 



Moreover, from equations (4.4) and (4.11), we can express the energy density (3.10) 
and pressure (3.11) of the effective stress energy tensor in the foreground space as 

b ( b b 2 \ . . 

p=- 3« 1 + 3« 2 - + « 3 4.14 

a V a or / 
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and 

b ( b\ ( b b 2 \b , 

p = — 2a 1 + a 2 - - «i + 2a 2 - + a 3 — T , (4.15) 
a \ a J \ a cr J a 

respectively. The g-space Bianchi-inspired conservation equation is then automatically 
satisfied for arbitrary b(t). 

Let us briefly consider the implications on the cosmological evolution of the g- 
space and /-space due to the inter-twining relation between the NECs. The simplest 
possibility is to have both NECs satisfied, which in this cosmological setting implies 
that the two cosmologies are conformally related, with the contribution to the effective 
stress-energy being an arbitrary cosmological constant that can be tuned at will. (It 
must be emphasized that, in principle, more general solutions are possible for particular 
combinations of the parameters appearing in L int .) 

Nevertheless, the case when the NECs are not saturated is more interesting, because 
either the foreground or background geometry has a NEC-violating effective stress- 
energy tensor. Assuming that the pathologies are relegated to the "other" gravitational 
sector, that is, in our own sector p > and p + p > 0, then, in view of equations (3.10), 
(3.13), (3.15) and (3.21), we have two possibilities: 

(i) Either p > and p + p < 0, which is phantom energy 

(ii) Or p < and p+p < 0, which corresponds to an even weirder kind of cosmological 
fluid known as dual dark energy [42]. 

Regarding (i), as is well known, phantom energy can easily lead to a big rip [43], big 
freeze [44, 45], future singularities [46], and other oddities [47-50]. The only possibility 
to avoid future doomsday in this scenario, under the assumption p + p > 0, seems to 
be by requiring p + p — > asymptotically (implying the same behaviour for p + p), 
leading to a cosmology which approximates a de Sitter behaviour at late times. In the 
second case, (ii), it must be noted that this effective dual dark energy is compatible 
with a Lorentzian geometry (3i/| = m 2 p/n + A — k/b 2 > 0), a positive cosmological 
constant A > (with A = Af + 0:4), or in a spatially hyperbolic universe, k < 0. As 
was studied in reference [51] for the case of a constant ratio w = p/p, a universe filled 
with dual dark energy and equipped with a positive cosmological constant (with k = 0) 
is completely regular and the A-term dominates at late times. Thus, in this case, the 
/-universe is asymptotically de Sitter at late times (for A > and any k). 

On the other hand, in massive gravity we are free to choose any arbitrary function 
b(t) since there are no equations of motion (2.10) for the background metric. Therefore, 
we can use the graviton mass terms to generate by hand an effective stress-energy 
contributing to cosmological fluid that has any properties we desire — the large scale 
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evolution of our universe is no longer governed by the matter in our universe, but 
can instead be forced to dance to the tune determined by the arbitrarily specified 
background geometry (see reference [19], and for earlier [pre-ghost-free] speculations 
along these lines see [3]). 

4.2 Massive-gravity black holes 

As we have already pointed out, the present ansatz can be motivated by generalizing the 
Gordon metric (and the relativistic acoustic metric) appearing in the field of analogue 
spacetimes [30, 31]. Following this spirit, and in view of the achievements of that field, 
one might expect this ansatz to be able to describe black hole geometries, at least 
when the metric describing the "medium" is the Minkowski metric. It must be noted 
that considering one of the metrics to be flat would implicitly restrict our attention to 
massive gravity, where one of the metrics is non-dynamical. This is the theory that we 
are implicitly assuming in the present subsection. 

Furthermore, let us note that a Minkowski medium in expression (3.3) implies 
9fiu = Vfiu, so 9fiu cannot be the dynamical metric, whereas it is Q^ u = f^ v that is 
dynamical. Usually in massive gravity the physical dynamical metric is denoted by g^ v 
and the background metric is f^, but this is simply a matter of notation. (Recall the 
symmetry of the interaction Lagrangian (2.8) which leads to the graviton contribution 
to the matter Lagrangian in massive gravity.) We maintain the notation of (3.3) for 
consistency with the analogue model programme. (In particular, all the results of this 
subsection would remain unchanged by adopting and specializing the "inverse ansatz" 
(3.4), considering g^ v to be the dynamical metric and f^ u = rj^ to be non-dynamical, 
and introducing an extra normalization factor in the four- velocity (4.19); since, in that 
case, it would be V that was a timelike unit vector with respect to the curved metric, 
while V would be a timelike unit vector with respect to the Minkowski metric.) 

It should be pointed out that assuming a flat background for massive gravity is not 
too strong an assumption when dealing with black holes, (though it would certainly 
be an overly stringent assumption for massive-gravity cosmology), at least as long as 
we are not presuming that both metrics can be written in a diagonal form in the same 
coordinate patch [28]. (See also [29].) In fact, to our knowledge, massive-gravity black 
holes (not bimetric gravity black holes) have been studied in the literature only by 
assuming flat backgrounds [22-25]. 

In massive gravity with a conformally flat background metric, our ansatz (3.1) can 
be written as 

= + (4.16) 

There are at least two ways of proceeding. 
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4.2.1 Isotropic coordinates 

Take V* = (1, 0), and set U 2 = ft 2 (l - £)■ Then 

ds 2 = -U 2 dt 2 + {l 2 \\dx\\ 2 . (4.17) 

This puts the metric into isotropic coordinates — and it is well known that this form 
of the metric is certainly sufficient (in principle) to deal with the Schwarzschild and 
Reissner-Nordstrom geometries, though the Kerr and Kerr-Newman geometries cannot 
be put in this form [52]. However in this coordinate system 7 = \fg~^f is singular (in 
the matrix sense) as any horizon (corresponding to U = 0) is approached. Thus once 
L int is "turned on" diagonal coordinate systems of this type behave unpleasantly - 
the stress tensor T tJ- u is then a singular matrix at the horizon, and other approaches 
might be more profitable. 

4.2.2 Horizon-penetrating coordinates 

Let us now adopt horizon-penetrating coordinates. Choose Q — 1, and express the 
background metric in the following gauge: 

ds 2 = rj^dx^dx" = -dt 2 + dr 2 + r 2 {d9 2 + sin 2 9 # 2 }, (4.18) 

Take the four-velocity to be 

^ = 7(1,-/3*); 7 M = - T (l,/3*); ^ = (/3 r , 0, 0). (4.19) 

Note that the usual relation between /3 and 7 is recovered when the normalization of 
V/j, with respect to 77^ is imposed. (The 7 appearing here should not be confused with 
the matrix 7 = \/g _1 f depending on the two metrics; its form would be that of the 
usual relativistic quantity. The intended meaning should be clear from context.) Thus, 
taking into account (4.16), (4.18) and (4.19), we can write 

ds) = - (1 - £ 7 2 ) dt 2 + 2^7 V7 2 - Idtdr + [1 + £ ( 7 2 - l)] dr 2 + r 2 dn 2 {2) . (4.20) 

Here dfl 2 2 ^ = d6 2 + sin 2 6 d(j) 2 must not be confused with the conformal factor Q. 
Let us now consider the specific function 

7 (r) = v 7 ! + K/r. (4.21) 

Then the foreground metric (4.20) can be written as 

ds 2 f = - (l - f - ^— ) dt 2 + 2 ^ Ki - K + r) dtdr + dr 2 + r 2 d{l 2 (4.22) 
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In its current form it is not canonically normalized at spatial infinity. If we now re-define 
the time variable as t = y/1 — 1 1, then 

ds 2 f = - ( 1 - , & ) dt 2 + 2 ^ K ^ r) dtdr +(l + ^)dr 2 + r 2 dn 2 (4.23) 
V (1 - r J y/l-Zr V r J 

which by explicit computation can easily be seen to be Ricci flat. This allows us to 
identify 

2m = (fT^ (424) 



and so 



2m\ - , J2m[2m{l-(i)+(ir}~ , / , 2m(l - 



ds=i = - ( 1 - — ] rft 2 + 2 v '"• L "" r — ^ ' ^' s dtdr+ I 1 + — ^ ) dr 2 +r 2 dtt 2 



(4.25) 

Now clearly in the limit £ — > I one recovers the Painleve-Gullstrand form of the 
Schwarzschild metric, whereas if £ = 0, then the black hole is expressed using the 
Eddington-Finkelstein metric. (For a related analysis, see [53].) 

One might think that there is some problem in performing a time re-definition 
once the gauge freedom is completely used to fix the form of the background metric 
(4.18). However, the same result can be obtained in a slightly different manner without 
explicitly performing any coordinate change. Let us start from the background metric 
already written using the new time coordinate. This is 

dsl = fj^dx^dx" = - 1 dP + dr 2 + r 2 {dd 2 + sin 2 9 d<p 2 }. (4.26) 

Considering now 

v» = 7 (VR, ; % = - 7 (^=|. > 

with 7 given by equation (4.21), and 

t = Vl#», (4.28) 

one can recover directly the foreground metric as given in (4.25). In this coordinate 
system the matrix 7 = -\/g~ l f and its inverse are non-singular (in the matrix sense) 
as any horizon is approached. Thus once L int is "turned on", the stress tensor T^ v is a 
non-singular matrix at the horizon, and can easily be tuned to be perturbatively small. 
Turning on the massive graviton should then not qualitatively disturb the horizon 
structure. 



(4.27) 
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Thus we see that the generalized Gordon ansatz is sufficiently flexible to deal with 
Schwarzschild black holes (and Reissner-Nordstrom black holes). This observation is 
important in that it gives us confidence that the ansatz is not too restrictive — the 
physics results we are obtaining above really are due to the interaction Lagrangian, not 
due to an overly restrictive metric ansatz. 

5 Generalized Kerr— Schild ansatz 

We now consider an ansatz which can be thought as being inspired by the Kerr-Schild 
metrics [54, 55]. This ansatz takes a similar form to that considered in the previous 
section, namely 



where the vector / M is now a null vector for foreground metric g^V'V = 0. So the 
geometries in question are conformally Kerr-Schild. It should be emphasized that the 
Kerr-Schild ansatz is an important and rather general one, which includes almost all 
spacetimes of physical interest [54, 55]. Note that (see figure 2): 

• For £ < 0, apart from a single null vector in common, the light cones of / lie 
strictly outside the light cones of g. 

• For £ = the light cones of / and g coincide; the two metrics are conformally 
related. (In this particular case the generalized Gordon and generalized Kerr- 
Schild ansatze coincide.) 

• For £ > 0, apart from a single null vector in common, the light cones of / lie 
strictly inside the light cones of g. 

For this ansatz, one has 



(5.1) 



(5.2) 



whence 




(5.3) 



and in general 




(5.4) 



The symmetric polynomials are given by 



ei( 7 ) = 4J2; 



e 2 ( 7 ) = 6ft 2 ; 



e 3 ( 7 ) = 4ft 3 . 



(5.5) 
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Figure 2. Light cones for the generalized Kerr-Schild ansatz: Depending on the value of the 
parameter £, and apart from one null ray in common, the foreground light cones lie strictly 
inside, on top of, or strictly outside the background light cones. 

(Note these quantities are independent of £.) Thus, taking into account equations 
(2.15) and (2.16), the stress-energy tensor can in this situation be written as 

T» v = ?l»l v + E8' l v , (5.6) 

where 

T = -^(ai + 2fia 2 + fi 2 « 3 ) , (5.7) 

and 

H = -0(3«! + 3fia 2 + a 3 tt 2 ). (5.8) 

Note that this stress-energy tensor can be interpreted as a null fluid (with null flux 
J 7 , and pressure H). This is a specialization of the "type II" stress-energy tensor 
considered by Hawking and Ellis [56]. Satisfaction of the NEC in this situation is 
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checked by contracting with a generic null vector (not necessarily P), and amounts to 
the condition 7 > 0. 

We can also write (2.15) 

r"„ = 7 1% + 3 (5.9) 

By inspection 

7 = 7; E = E + L int = n( ai + 3a 2 n + a 3 n 2 ). (5.10) 

Finally, 

T u = 7l% + Z8» u . (5.11) 
Noting that now y/—f = Q 4 y/—g, we have 

7=-^P, S=-ijg. (5.12) 
Regarding the NEC the key observation is 

T— (5.13) 
We again very explicitly see how foreground and background NECs are anti-correlated [14]. 

6 Discussion 

In the technical analysis of ghost-free bimetric gravity a central role is played by the 
matrix square root 7 = \A? -1 / relating the foreground and background metrics. Mo- 
tivated by the fact that this object can be somewhat awkward to deal with, we have 
developed two lines of attack. In a companion paper we have developed a number of 
mathematical results that help us manipulate this quantity [14], ultimately leading to a 
rather general analysis of the interplay between the foreground and background NECs. 
In the current article we have taken a slightly different tack, and have characterized the 
circumstances under which it is possible to simply and explicitly compute this matrix 
square root. 

In order to simplify the problem we have assumed that the metrics are related in 
a very particular form, fulfilling what we called the "generalized Gordon ansatz" . This 
ansatz can be motivated from the analogue spacetime programme, and has been proven 
to be amazingly powerful in this paper. It implies that the effective stress-energy tensor 
automatically takes the form of a perfect fluid. This ansatz is extremely natural in a 
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cosmological setting, and helps one understand why bimetric FLRW cosmologies have 
proved to be relatively tractable. 

Moreover, it has also allowed us to recover easily the relation between the NECs 
in the two spaces [14] . That is, we have shown that the simultaneous fulfillment of the 
NEC associated to the modifications of general relativity in both foreground and back- 
ground spaces is only possible if it saturates; that is if the effect of these modifications 
are equivalent to the consideration of a cosmological constant. 

In counterpoint to these rather strong conclusions, it should be noted that the 
generalized Gordon ansatz seems sufficiently general to study the most common physical 
situations one might be interested in. In fact, we have shown that it can describe the 
general class of cosmological solutions in bimetric gravity, as well as black holes in 
massive gravity. 

We have also considered a less general ansatz based on the Kerr-Schild metrics. 
The effective stress-energy tensor in this case has not quite so natural a form; it is 
equivalent to a null fluid. Nevertheless, even in this case similar conclusions can be 
extracted about the simultaneous fulfillment of both NECs — this is only possible in 
the trivial case in which both metrics have the same light cone structure. 
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